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ABSTRACT 


The solar air heater has an important place among solar heat collectors. It can 
be used as sub-systems in many systems meant for the utilization of solar 
energy. Possible applications of solar air heaters are drying or curing of 
agricultural products, space heating for comfort regeneration of 
dehumidifying agents, seasoning of timber, curing of industrial products 
such as plastics. When air at high temperature is required the design of a 
heater becomes complicated and very costly. As far as the ultimate application 
for heating air to maintain a comfortable environment is concerned, the solar 
air heater is the most logical choice. In general solar heaters are quite suitable 
for low and moderate temperatures application as their design is simple. 
These solar air heaters have low heat transfer efficiency which can be 
improved using geometrical modifications like optimizing duct geometry or 
adding artificial roughness. The use of artificial roughness on the underside of 
the absorber plate is an effective and economic way to improve the thermal 
performance of a solar air heater. Several experimental investigations, 
involving different types of roughness elements, have been carried out to 
improve the heat transfer from the absorber plate to air flowing in solar air 
heaters. In this paper the CFD analysis on heat transfer and friction in 
rectangular ducts roughened with broken 'S' shaped ribs has been presented. 
The relative gap width [g/e] is varied from 0.5 to 2.5 and other parameter are 
constant. The effects of gap width [g/e] on Nusselt number, friction factor 
and thermo-hydraulic performance parameter have been discussed and 
results compared with smooth duct under similar conditions. It is found that 
the maximum heat transfer and friction characteristic at a relative gap width 
of 1. 
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I. 


INTRODUCTION 


In general solar heaters are quite suitable for low and 
moderate temperatures application as their design is simple 
technique. It has been found that the main thermal 
resistance to the convective heat transfer is due to the 
formation of boundary layer on the heat transferring surface. 
Efforts for increasing heat transfer have been directed 
towards artificial roughness on a surface is an effective 
technique to enhance heat transfer to fluid flowing in the 
duct by Bhagoria JLet al.[l].Artificial roughness in the form 
of wires and in various arrangement has been used to create 
turbulence near the wall or to break the boundary layer. 

Thus the artificial roughness can be employed for the 
enhancement of heat transfer coefficient between the 
absorber plate and air, improving the thermal performance 
of solar air heater contrary to this, the excessive disturbance 
to the boundary layer creates more friction resulting in more 
pumping power by Bhagoria JLet al [2] .Application of 
artificial roughness in solar air heater owes its origin to 
several investigations carried out for enhancement of cooling 
of turbine blades passage. Several investigations have been 


carried out to study effect of artificial roughness on 
heat transfer and friction factor for two opposite roughened 
surfaces by Han[3, 4], Han et al.[5-7], Wrieght et al.[8], 
Lueetal, [9-11], Taslim et al. and Hwang[14], Han and 
Park[15], Park et al.[16] developed by different 
investigators. The orthogonal ribs i.e. ribs arranged normal to 
the flow were first used in solar air heater and resulted in 
better heat transfer in comparison to that in conventional 
solar air heater by Prasad k, Mullick S. C. et al [17].Many 
investigators Gao x sunden B[18], Han J.C, Glicksman LR, 
Rohsenow WM[19], Prasad BN, Saini JS{20], Taslim ME, 
Li T, Kercher Dm[21], Webb RL, Eckert Erg, Goldstein 
RJ[22] have reported in detail the Nu and f for orthogonal and 
inclined rib-roughened ducts. The concept of V-shaped ribs 
evolved from the fact that the inclined ribs produce 
longitudinal vortex and hence higher heat transfer. In 
principal, high heat transfer coefficient region can be 
increased two folds with V-shape ribs and hence result in 
even higher heat transfer et al. [21].The beneficial effect on 
Nu and f caused by V-shaping of ribs in comparison to angled 
ribs has been experimentally endorsed by several 
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investigators Geo X, Sunden B{23], Karwa R.[24], Kukreja 
RT, LueSC, McMillin RD[25], LauSC, McMillin RD, Han 
JC[26], for different roughness parameters and duct aspect 
ratios. For V-shape ribs, the inter-rib local heat transfer 
coefficient reduces from leading edge(s] to trailing edge(s] in 
transverse direction[21, 23, 24], However in the flow 

direction, the inter- rib local heat transfer coefficient varies 
like saw tooth [21, 23, 24].In addition, multiple V-ribs 
have also been investigated with the anticipation that the 
more number of secondary flow cells may result in still 
higher heat transfer et at Lanjewar A, Bhagoria JL, Sarviya 
RM[27], Hans VS, SainiRP, Saini JS[28]. The second law of 
thermodynamic analysis combined with a standard design 
procedure of a thermal system gives us invaluable insight 
into the operation of the system. The studies of Han et al.[5], 
Lau et al.[10] and Taslim et al.[12] not covered the wide 
range of roughness and operating parameters as would be 
required for detailed analysis for detailed optimal design or 
selection of roughness parameter to be used in conventional 
solar air heaters. Most of the investigations carried out so far 
have applied artificial roughness on two opposite wall with 
all four walls being heated. However in case of solar air 
heater, roughness elements are applied to heated wall while 
remaining three walls are insulated. Heated wall consists of 
absorber plate and is subjected to uniform heat flux 
(insolationJ.This makes fluid flow and heat transfer 
characteristics distinctly different from those found in case of 
two roughened wall and four heated wall duct. 

Heated wall consists of absorber plate and is subjected to 
uniform heat flux(insulation].This makes fluid flow and heat- 


transfer characteristics distinctly different from those found 
in case of two roughened walls and four heated wall duct 
Producing a gap in the inclined rib is found to enhance the 
heat transfer by breaking the secondary flow and producing 
higher level of turbulence in the fluid downstream of the rib 
.A similar gap in both the limbs of v-rib further enhances the 
heat transfer by introducing similar effects in both the limbs. 
Further the use of multi v-rib across the width of the plate is 
found to enhance the heat transfer by increasing the number 
of secondary flow cells several times. It is thought that 
producing gaps in all the limbs of multi-v geometry will 
bring about considerably large enhancement in comparison 
to that of simple single v-rib arrangement. It will therefore 
be pertinent to investigate the effect of various geometrical 
and flow parameters on the heat transfer and friction 
characteristics of rectangular duct having its absorber plate 
roughened with 'S 'shaped with gap. 

II. Computational Fluid Dynamics 

Computational fluid dynamics or CFD is the analysis of 
systems involving fluid flow, heat transfer and associated 
phenomena such as chemical reactions by means of 
computer-based simulation. The technique is very powerful 
and spans a wide range of industrial and non-industrial 
application areas. The 2-dimensional solution domain used 
for CFD analysis has been generated in ANSYS version 14.5 
(workbench mode] as shown in Fig.l.The solution domain is 
a horizontal duct with broken 'S' shaped ribs roughness on 
the absorber plate at the underside of the top of the duct 
while other sides are considered as smooth surfaces. 



Fig.l. Showing the geometric dimension of the working model 



Complete duct geometry is divided into three sections, 
namely, entrance section, test section and exit section. A 
short entrance length is chosen because for a roughened 


duct, the thermally fully developed flow is established in a 
short length 2-3 times of hydraulic diameter. The exit 
section is used after the test section in order to reduce the 


@ IJTSRD I Unique Paper ID - IJTSRD33610 | Volume - 4 | Issue - 6 | September-October 2020 Page 1338 






























International Journal of Trend in Scientific Research and Development [IJTSRD] @ www.ijtsrd.com elSSN: 2456-6470 


end effect in the test section. The top wall consists of a 0.5 
mm thick absorber plate made up of aluminum. Artificial 
roughness in the form of small diameter galvanized iron (G.I] 
wires is considered at the underside of the top of the duct on 
the absorber plate to have roughened surface, running 
perpendicular to the flow direction while other sides are 
considered as smooth surfaces. A uniform heat flux of 1000 
w/m2 is considered for computational analysis. 



Fig.3 Meshing of computational Domain for broken 'S' 
ribs Roughness 


A non-uniform mesh is shown in Fig.SPresent mesh 
contained 191, 061 quad cells with non-uniform quad grid 
of 0.21 mm cell size. This size is suitable to resolve the 
laminar sub-layer. For grid independence test, the number 
of cells is varied from 113, 431 to 207, 147 in five steps. It 
is found that after 191, 0481 cells, further increase in cells 
has less than 1% variation in Nusselt number and friction 
factor value which is taken as criterion for grid 
independence. 

In the present simulation governing equations of continuity, 
momentum and energy are solved by the finite volume 
method in the steady-state regime. The numerical method 
used in this study is a segregated solution algorithm with a 
finite volume-based technique. The governing equations are 
solved using the commercial CFD code, ANSYS Fluent 14.5. 
No-slip conditions for velocity in solid surfaces are assumed 
and the turbulence kinetic energy is set to zero on all solid 
walls. The top wall boundary condition is selected as constant 
heat flux of 1000 W/m^ and bottom wall is assumed at 
adiabatic condition. A uniform air velocity is introduced at 
the inlet while a pressure outlet condition is applied at the 
exit. The Reynolds number varies from 2000 to 16000 at the 
inlet. The mean inlet velocity of the flow is calculated using 
Reynolds number. Constant velocity of air is assumed in the 
flow direction. The temperature of air inside the duct is also 
taken as 300 K at the beginning. At the exit, a pressure 
outlet boundary condition is specified with a fixed pressure 
of 1.013x105 Pa. 

III. RESULTS AND DISCUSSION 
A. Heat Transfer Cbaracteristics and Friction Factor 
Characteristics 

Fig.4 shows the effect of Reynolds number on average 
Nusselt number for different values of relative gap width 
(g/e] and fixed other parameter. The average Nusselt 
number is observed to increase with increase of Reynolds 
number due to the increase in turbulence intensity caused by 
increase in turbulence kinetic energy and turbulence 
dissipation rate. 
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Fig. 4 Variation of Nusselt number with Reynolds 
number for different Values of relative gap width (g/e) 

Effect of the relative gap width (g/e] on heat transfer is also 
shown typically in Fig. 4. It can be seen that the enhancement 
in heat transfer of the roughened duct with respect to the 
smooth duct also increases with an increase in Reynolds 
number. It can also be seen that Nusselt number values 
increases with the increase in relative gap width (g/e] of up 
to 1 and then decrease for a fixed value of roughness pitch 
(P]. The roughened duct having broken 'S' shaped with 
relative gap width (g/e] of 1 provides the highest Nusselt 
number at a Reynolds number of 16000. For circular rib the 
maximum enhancement of average Nusselt number is found 
to be 2.45 times that of smooth duct for relative gap width 
(g/e] of 1.0 at a Reynolds number of 16000. 

The heat transfer phenomenon can be observed and 
described by the contour plot of turbulence intensity. The 
contour plot of turbulence intensity for broken 'S' shaped ribs 
is shown in Fig.5 (a, b and c].The intensities of turbulence 
are reduced at the flow field near the rib and wall and a high 
turbulence intensity region is found between the adjacent 
ribs close to the main flow which yields the strong influence 
of turbulence intensity on heat transfer enhancement. 
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The value of this parameter higher than unity ensures that it 
is advantageous to use the roughened duct in comparison to 
smooth duct. The thermo-hydraulic parameter is also used to 
compare the performance of number of roughness 
arrangements to decide the best among these. The variation 
of thermo-hydraulic parameter as a function of Reynolds 
number for different values of relative gap width (g/e] and 
investigated in this work has been shown in Fig. 7. For all 
values of relative gap widths, value of performance 
parameter is more than unity. Hence the performance of solar 



Figure. 6 also shows that the friction factor decreases with 
the increasing values of the Reynolds number in all cases as 
expected because of the suppression of laminar sub-layer for 
fully developed turbulent flow in the duct. It can also be seen 
that friction factor values increase with the increase in 
relative gap width (g/e] up to 1 and then decrease for fixed 
value of roughness pitch, attributed to more interruptions 
in the flow path. 

B. Thermo-Hydraulic Performance 

It has also been observed from Figures 4 and 6 that the 
maximum values of Nusselt number and friction factor 
correspond to relative gap width of 1.0, thereby, meaning 
that an enhancement in heat transfer is accompanied by 
friction power penalty due to a corresponding increase in the 
friction factor. Therefore, it is essential to determine the 
effectiveness and usefulness of the roughness geometry in 
context of heat transfer enhancement and accompanied 
increased pumping losses. In order to achieve this objective, 
Webb and Eckert [23] proposed a thermo-hydraulic 
performance parameter 'p', which evaluates the 

enhancement in heat transfer of a roughened duct compared 
to that of the smooth duct for the same pumping power 
requirement and is defined as. 


Thermal enhancement factor = ■ 


Nu/ 

/] 


Nu, 




and relative roughness pitch on the heat transfer coefficient 
and friction factor have been studied. In order to validate the 
present numerical model, results have been compared with 
available experimental results under similar flow conditions. 
CFD Investigation has been carried out in medium Reynolds 
number flow (Re = 2000-16, 000]. The following 

conclusions are drawn from present analysis: 

1. The Renormalization-group (RNG] k-e turbulence model 
predicted very close results to the experimental results, 
which yields confidence in the predictions done by CFD 
analysis in the present study. RNG k-e turbulence model 
has been validated for smooth duct and grid 
independence test has also been conducted to check the 
variation with increasing number of cells. 

2. The roughened duct having broken 'S' shaped rib with 
relative gap width of 1.0 provides the highest Nusselt 
number at a Reynolds number of 16000. 

3. For rectangular rib the maximum enhancement of 
average Nusselt number is found to be 2.46 times that of 
smooth duct for relative gap width of 1.0 at a Reynolds 
number of 11000. 

4. The roughened duct having broken 'S' shaped rib with 
relative gap width of 1.0 provides the highest friction 
factor at a Reynolds number of 3500. 

5. For broken 'S' shaped rib the maximum enhancement of 
average friction factor is found to be 3.48 times that of 
smooth duct for relative gap width of 1.0. 


@ IJTSRD I Unique Paper ID - IJTSRD33610 | Volume - 4 | Issue - 6 | September-October 2020 Page 1340 
































International Journal of Trend in Scientific Research and Development [IJTSRD] @ www.ijtsrd.com elSSN: 2456-6470 


It is found that the thermal hydraulic performance of relative 

gap width of 1.0 is maximum. 
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